Topological defects: A problem for cyclic universes? 
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We study the behaviour of cosmic string networks in contracting universes, and discuss some of 
their possible consequences. We note that there is a fundamental time asymmetry between defect 
network evolution for an expanding universe and a contracting universe. A string network with 
negligible loop production and small-scale structure will asymptotically behave during the collapse 
phase as a radiation fluid. In realistic networks these two effects are important, making this solution 
only approximate. We derive new scaling solutions describing this effect, and test them against 
high-resolution numerical simulations. A string network in a contracting universe, together with 
the gravitational radiation background it has generated, can significantly affect the dynamics of the 
universe both locally and globally. The network can be an important source of radiation, entropy 
and inhomogeneity. We discuss the possible implications of these findings for bouncing and cyclic 
cosmological models. 



I. INTRODUCTION 

Cosmological scenarios involving oscillating or cyclic 
universes have been know for a long time with in- 
terest in them varying according to the latest theoretical 
prejudices or observational constraints. Recent interest 
has been associated with a cyclic extension of the ekpy- 
rotic scenario 0]- A related result was the realization 
1^ that the presence of a scalar field seems to be 
necessary to make cosmological scenarios with a bounce 
observationally realistic. And if scalar fields are present, 
then one should contemplate the possibility of topological 
defects being formed. 

It is thought that the early universe underwent a se- 
ries of phase transitions, each one spontaneously break- 
ing some symmetry in particle physics and giving rise to 
topological defects of some kind 0, Q, which in many 
cases can persist throughout the subsequent evolution of 
the universe. 

In the present work we study cosmic string evolution 
in a collapsing universe, following up on and generaliz- 
ing the results of 0, and discuss in much greater detail 
some implications of the presence of cosmic strings (and 
cosmic defects in general) for bouncing universes. In a 
bouncing universe scenario the properties of the universe 
in the expanding phase depend on physics happening in 
a previous collapsing phase (before the bounce). For 
this reason, if defects do exist in these models, it is cru- 
cial to understand their evolution and consequences in 
both the expanding and collapsing phases. Up to now 
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all these studies, be they analytic 0, 0, 0, 0, or 
numerical fl^L H^ . have only been undertaken for 
the expanding case, and it is clear that while some re- 
sults may be expected to carry over to the contracting 
phase, some others clearly won't. This will have con- 
sequences not only for the standard string seeded (or 
hybrid) structure formation scenario but also for 
other 'non-standard' scenarios involving defects, such as 
the production of adiabatic and nearly Gaussian density 
fluctuations [iM , or those involving anisotropic or in- 
homogeneous universes |^ . 

In particular, we expect that cosmic strings will be- 
come ultra-relativistic, behaving approximately like a ra- 
diation fluid. This means that a cosmic string network, 
both directly and through the gravitational radiation 
emitted by the small loops it produces, will soon become 
a significant source of entropy (and also of inhomogene- 
ity) . Hence a cosmic string network is a further problem 
for cyclic universes if a suitable and efficient mechanism 
for diluting the entropy is not available. 

We should point out at the outset that if/when dur- 
ing the collapse phase one reaches the Hagedorn tem- 
perature, one expects the string network to quickly dis- 
solve. However this is largely irrelevant for the points 
being made in this paper: the radiation, entropy and 
anisotropy produced by the network will obviously still 
be left behind if the network does dissipate. On the other 
hand, it need not be the case that the Hagedorn temper- 
ature is reached or, more specifically, that the collapse 
has to continue all the way to the 'Big Crunch' (where 
there is no known sensible description of the physics in- 
volved) . Cosmological models do exist where the bounce 
takes place at finite size. Indeed, such models seem to 
be relatively common in scenarios with extra dimensions 
and scalar fields, although this issue is somewhat debat- 
able. 

The outline of this paper is as follows. In Sect. UTI we 
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present a very brief overview of previous work on cyclic 
universes, and discuss further motivation for our work. 
Then in Sect lllll we introduce the basic dynamical prop- 
erties of cosmic string networks, and after a 'warm-up' 
example of a circular string loop we successively discuss 
analytic scaling solutions in various regimes in a contract- 
ing universe. In Sect. IIVI we describe high-resolution 
numerical simulations of string networks in contracting 
universes using the AUen-Shellard string code ^5|, and 
then compare and contrast our analytic and numerical 
results. Possible cosmological consequences of these re- 
sults are discussed in Sect. and finally we present 
some conclusions in Sect. I VII Throughout the paper we 
will use units in which c = h = 1. 



II. AN OVERVIEW OF BOUNCING 
COSMOLOGICAL MODELS 

Oscillating universes arise naturally as classical exact 
solutions of the Einstein equations, and have been ex- 
plored in several contexts in an attempt to solve some 
of the cosmological enigmas. For example, they could 
conceivably solve the flatness and horizon problems. For 
oscillating universes whose size at maximum expansion 
increases in each cycle the flatness problem may be solved 
because it is ever less likely for an observer to find himself 
in a non-flat region. Provided there is a causal correla- 
tion for the microphysics at the bounce, the age for cyclic 
universes may also be large enough to solve the horizon 
problem. 

However, oscillating models are not without their 
caveats. The first of these was pointed out a long time 
ago by Tolman P|. By assuming a closed universe with 
zero cosmological constant in the context of the General 
Relativity theory, he noted that entropy is generated at 
each cycle and so the total entropy in the universe grows 
from cycle to cycle. Hence the period of each cycle is 
larger that the previous one, and extrapolating back in 
time one finds that the sum is finite. In other words, 
the universe must still have had a beginning. These 
points were also discussed by Rees , and formalized by 
Zel'dovich and Novikov 23] and subsequently by others 

MM 

A thorough analysis of oscillating universe solutions 
was first carried out by Barrow and Dabrowski 26] . They 
showed that in closed universes filled with a perfect fluid 
of matter and/or radiation, with total entropy increas- 
ing from cycle to cycle, a positive cosmological constant 
or any non-zero stress violating the strong energy con- 
dition will eventually halt the oscillations after a finite 
number of cycles. In particular when the cosmological 
constant dominates the universe approaches the de Sit- 
ter solution. Dabrowski subsequently generalized the 
analysis to include negative pressure matter which vio- 
lated the strong energy condition, and also provided a 
somewhat simplistic description of the scalar field case. 

Looking for an alternative model to inflation that could 



explain the scale invariance of the observed power spec- 
trum of fluctuations, Durrer and Laukenmann pre- 
sented a model of a closed universe dominated by radi- 
ation or with an intermediate matter-dominated period 
where small black holes could have been formed. They 
postulated that density perturbations produced gravita- 
tional entropy which would have to be transformed into 
radiation entropy at the bounce. This could explain most 
of the present day radiation entropy without having over- 
production, provided one also assumed that perturba- 
tions never became strongly non-linear in the previous 
cycle (as the mass fluctuation would in fact generally di- 
verge at the big crunch). 

More recently, a detailed analysis by Peter and Pinto- 
Neto has shown quite generically that no detectable 
bounce is observationally allowed in any universe which 
around the epoch of the bounce is described by Einstein 
gravity together with hydrodynamical fluids, since un- 
der those conditions scalar density perturbations would 
become non-linear well before nucleosynthesis. However, 
this 'no-bounce conjecture' can in some circumstances be 
evaded by adding a free scalar field which can dominate 
the universe during the bounce . 

One of the earliest discussions of bouncing universes 
with scalar fields was by Hawking in the context of the 
no-boundary proposal |23; of course, many of the dis- 
cussions noted earlier are derivative of Hawking's cosmo- 
logical singularity theorems. Early in the development 
of quantum cosmology, he studied a simple chaotic in- 
flation model and suggested that his boundary condi- 
tions selected only a special set of time-symmetric and 
eternally bouncing trajectories. Subsequent work, how- 
ever, showed that these special bouncing states were not 
generic [s^, IsJ], with typical universes diverging from 
what are unstable trajectories [s^] and re-collapsing to 
singularities. 

Kanekar et al. Q have also considered a massive scalar 
field in a closed universe for the Einstein theory. This 
is a possible alternative mechanism which can still lead 
to an increase in the volume of an oscillating universe 
without requiring entropy production. They suggest that 
the oscillations of the universe force the scalar field to 
also oscillate about the minimum of its potential. The 
work done by or on the scalar field creates an asymmetry 
between the expansion and collapse epochs which can 
result in the increasing volume of the universe. It was 
also argued that the presence of other matter fields will 
not change this picture, provided that the interactions 
between the matter and scalar fields are sufficiently weak. 

Finally, we should mention the cyclic extension 2] of 
the ekpyrotic model |U |3 . Here the contraction and 
expansion phases correspond to the epochs before and af- 
ter the collision of two thrcc-branes along a fifth dimen- 
sion. During the expansion there are successive periods 
of radiation, matter and vacuum domination (the latter 
of which can dilute the entropy, black holes and other 
debris produced during each cycle), and the inter-brane 
potential can be carefully chosen to ensure the cyclic be- 
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haviour. Hence the model can potentially solve many 
of the cosmological issues raised above. However, it ap- 
pears questionable whether density perturbations in this 
model, matched across the bounce, can reproduce the 
Harrison-Zel'dovich spectrum, as has been pointed out 
by a number of authors [sE HE 1^ 1^ • 

In the present work, motivated by the Kanekar et al. 
results, we replace the scalar field by a cosmic string 
network, which is also expected to display an asym- 
metric behaviour between the contraction and expansion 
epochs. In particular, while during expansion a cosmic 
string network will quickly evolve towards a linear scal- 
ing regime ^E (except in very particular circumstances) , 
we shall see that this is not the case during a phase of 
collapse. Indeed, in this case, a string network where 
loop production and small-scale 'wiggles' were negligible 
would asymptotically behave like a radiation fluid. As 
it turns out, these two effects are relevant throughout 
the collapse (in fact, even more so than in the expanding 
phase), and this solution is only approximate. 

On the other hand, there is further radiation being 
produced: since cosmic string loops decay gravitation- 
ally, the enhanced loop production will give rise to an 
enhanced gravitational radiation background. All in all, 
a cosmic string network will add a significant contribu- 
tion, in the form of radiation, to the energy (and hence 
also entropy) budget of a contracting universe, which will 
become ever more important as the contraction proceeds. 
In the following sections we will demonstrate this be- 
haviour and discuss its cosmological consequences. 



III. COSMIC STRING EVOLUTION: BASICS 
AND ANALYTIC METHODS 

In this section we discuss the basic dynamical proper- 
ties of cosmic string networks, using both analytic and 
numerical methods, and use these tools to characterize 
their behaviour in bouncing universes. Both the ana- 
Mjc^lM 111 S EE Ef El El and the numerical 
[lEl lid Il7j| tools which we shall use rely on previous 
work by some of the present authors. In what follows 
we will limit ourselves to describing the features that are 
directly relevant for our analysis. We refer the reader to 
the original references for a more detailed discussion and 
derivations of some key results. 



A. Basics of string dynamics 

In the limit where the curvature radius of a cosmic 
string is much larger than its thickness, we can describe 
it as a one-dimensional object so that its world history 
can be represented by a two-dimensional surface in space- 
time (the string world-sheet) 

x" = {a") ; a = 0,1; z/ = 0,l,2,3 (1) 



obeying the usual Goto-Nambu action. For simplicity we 
start by considering the string dynamics in a flat FRW 
universe with line element 

ds^ ^a^(T])(df]^ ~dx^). (2) 

Identifying conformal and world-sheet times and impos- 
ing that the string velocity be orthogonal to the string 
direction (i.e. x • x' = 0) the string equations of motion 
take the form 

x + 27^ (l-i^) i = 1 J^^^ (3) 

e + 2n±^e = (4) 

where the 'coordinate energy per unit length', e, is de- 
fined by 



Ti. — d/a, and dots and primes are derivatives with re- 
spect to the time and space coordinates. Note that if the 
universe is not spatially flat then the microscopic string 
equations of motion 13141) will have additional curvature 
corrections — see |121I131I42| . 

The evolution of the scale factor is described by the 
Friedmann equation 

+ K = nlinl,a-^ + n°a-^ + nla^) (6) 

where we have considered only three contributions to the 
energy density of the universe, namely matter, radiation 
and a cosmological constant. Note that this assumes that 
the cosmic string density is comfortably subdominant 
and need not be included in the Friedmann equation. 
Although this is usually the case in most cosmological 
scenarios, it need not happen all the time. 

B. A simple example: the circular loop 

We will now consider the evolution of a circular cosmic 
string loop in a cyclic universe. Although realistic loops 
chopped off by the network are of course highly irregular 
(and, crucially, have a considerable amount of small-scale 
'wiggles' 0,[lEE3l): the circular solution is still a useful 
starting point and can provide some intuition for things 
to come |39j . In this example we will assume for simplic- 
ity that the universe is spatially flat and that there is a 
negative cosmological constant which makes the universe 
collapse. 

The loop's trajectory can be described by 

X = r{?]) (sin 9, cos 6, 0) . (7) 

with < 9 < 27T. We define R = \r\aj as the 'invariant' 
loop radius (which is proportional to the energy of the 
loop), and 7 = (1 — ?;^)~^/^, with v = \r\ being the 
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microscopic speed of the loop. The microscopic string 
equations of motion (|3I4|I then become 



r = (l-r^) [ - 2Hr 



(8) 



where H — Ti./a, with the velocity equation being obvi- 
ously obtainable form this, or alternatively in terms of 
the invariant loop radius 



dR 
'dt 



(1 - 2v^)HR. 



(9) 



This latter equation coincides with the averaged evo- 
lution equation for R, obtained in the context of the 
velocity-dependent one-scale model j3 [sp , while for 
the averaged velocity equation this model yields 



dv 
'di 



(1 



k{v) 



2Hv 



(10) 



where k(v) is the momentum parameter which is thor- 
oughly discussed in [T3|. 

Simple analytic arguments show that a loop whose ini- 
tial radius is much smaller than the Hubble radius will 
oscillate freely with a constant invariant loop radius and 
an average velocity v ~ l/\/2. (Note that we are assum- 
ing units in which c = 1.) On the other hand, once the 
collapse phase begins, we will eventually get to a stage 
in which the physical loop radius becomes comparable to 
the Hubble radius ar ~ and then gets above it. In 
this regime the loop velocity is typically driven towards 
unity V —^ 1 and it is straightforward to show that the 
invariant loop length grows as i? oc and the Lorentz 
factor as 7 oc a~^. Despite its growing energy R, the 
actual physical loop radius ar = R/"f oc a, so the loop 
shrinks with the scale factor and inexorably follows the 
collapse into final big crunch singularity. 

Importantly, note that this relativistic final state for 
a loop in a collapsing universe is generic and quite dif- 
ferent to the initial condition usually assumed for super- 
horizon loops in the expanding phase (created, say, at 
a phase transition). In the latter case, the loops begin 
with a vanishing velocity which only becomes significant 
when each of them falls below the Hubble radius. Such 
evolution cannot be reproduced in reverse during the col- 
lapsing phase without fine-tuning the velocity as the loop 
crosses outside the Hubble radius. This simple fact intro- 
duces a fundamental time asymmetry for string evolution 
in a cyclic universe. In fact, something analogous hap- 
pens for all other topological defects. 

The analytic expectations for our circular loop solution 
have been confirmed by a numerical study, see Fig. ^ 
Here we plot both the microscopic and averaged loop 
sizes and velocities, for four loops of different sizes in 
a universe filled with matter (the radiation case being 
analogous, except for a slight difference in the timescale 
required to reach this asymptotic regime). Obviously a 
loop which has a considerable size relative to the horizon 
to begin with will not oscillate. We can see that the 
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FIG. 1: Comparing the microscopic (darker curves) and av- 
eraged (lighter curves) evolution for four circular loops in a 
matter-dominated, closed universe. The conformal time was 
chosen so that the big bang occurs at ?7 = and the big crunch 
at r; = 1. Displayed are the loop velocity v (top left) and in- 
variant radius R (top right), as well as the ratio of the loop 
and Hubble radii R/dn (bottom left) and the product of the 
loop radius and the scale factor Ra (bottom right), showing 
the asymptotic R oc behaviour. 



averaged quantities provide a very good description of 
the dynamics. 

The main caveat with this solution is that it can not 
account for the presence of small-scale structures, whose 
build up we in fact expect to be enhanced in the collaps- 
ing phase. One expects that the presence of small-scale 
wiggles can significantly delay the onset of this asymp- 
totic regime — an expectation that we will confirm in what 
follows. 



C. String network evolution: analytic expectations 

Two different but complementary approaches are avail- 
able to study the evolution of a cosmic string network: 
one can resort to large numerical simulations [3 Q 
(which are intrinsically difficult and time consuming, as 
one is dealing with highly non- linear objects), or one can 
develop analytic tools 9^, JT, 12, JJi, J£, ^42,] which pro- 
vide an averaged (or 'thermodynamical') description of 
the basic properties of the network. In what follows we 
shall briefly describe the best motivated of these analytic 
models, the velocity-dependent one-scale (VOS) model 
[Hill 111 m, and try to use it to deduce the basic 
properties of a cosmic string network during a phase of 
contraction. In the following subsection we will test these 
against numerical simulations. 
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1. The VOS model 



The VOS model describes the string dynamics in terms 
of two 'thermodynamical' parameters: the string RMS 
velocity, v^o, defined by 



9 9 [x.^eda 

vL =< >= 



Jeda 



(11) 



and a single length scale, L, which can be variously inter- 
preted as the long string correlation length, its curvature 
radius, or simply a measure of the energy density (see 
below). It is important to note that in the context of 
this model, all these three length scales are assumed to 
be identical • Of course this assumption is not always 
realistic and models exist where it is relaxed 10, 12, 43| . 
but it must be assumed if one wants to use the model 
in this formulation. The string network is thus assumed 
to be a Brownian random walk on large enough scales, 
characterized by a correlation length L. Hence one can 
simply relate it with the energy density in long strings as 



Po 



L2 ' 



(12) 



where ii = fiQ is the string mass per unit length. Note 
that given that we will be considering relativistic veloci- 
ties, the commonly used 'correlation length' L is really a 
measure of the invariant string length or energy, rather 
than the typical curvature radius of the strings. By in- 
cluding the appropriate Lorentz factor 700 — (1— i;^)~^/^ 
for the long strings, we can denote the physical correla- 
tion length by 



(13) 



Note also that in defining (|12|) we are taking the string 
mass per unit length to be a constant. More generally, 
if one wanted to explicitly account for the presence of 
small-scale 'wiggles' one would need to define a varying 
'renormalized' mass per unit length 'T^ |43| . which would 
then be related to the 'bare' mass by 



(14) 



With this assumption the VOS model has one phe- 
nomenological parameter c, commonly called the loop 
chopping efficiency, which describes the rate of energy 
transfer from the long-string network to loops. The evo- 
lution equations then take the following form [iTl Is^ 



+ 8TG^lvl, (15) 



dvr, 



dt 



L 



- 2Hv„ 



(16) 



The final term in the evolution equation for the corre- 
lation length describes the effect of gravitational back- 
reaction. Notice that we are not including in either equa- 
tion additional terms arising from friction due to parti- 
cle scattering 13j, which could conceivably be important 



FIG. 2: The momentum parameter k{v) in the velocity- 
dependent one-scale string evolution model. Notice that it 
vanishes at v'^ — 1/2 and v'^ — 1, and it is negative in be- 
tween these two values. 



during the final stages of collapse. We shall return to 
this point below. Here k{voo) is the momentum param- 
eter, which is thoroughly discussed in 0, and whose 
dependence on the string velocity. 



k{v) 



2V2 



2\/2z;3)l 



1 + 



(17) 



is shown in Fig. |21 Notice that this is positive for < 
< 1/2 and negative for 1/2 < v"^ < 1; this turns out 
to be crucial for some of what follows. 

As a final remark, we point out that these evolution 
equations are approximate, and are neglecting higher- 
order terms |12|. It is in principle possible that in this 
particular case ultra-relativistic corrections may become 
significant, for example there could be powers of 7 in 
various terms above. On the other hand, one physically 
expects that the velocity increase will also be accompa- 
nied by an enhanced build-up of small-scale 'wiggles', for 
which one might need further corrections [T^ |43 | . ^From 
a 'phenomenological' modelling point of view, the latter 
are much harder to calculate than the former, and as far 
as one knows this can only be done in the context of 
proper wiggly models of string evolution [i^l l43 | . Thus 
in the present work we chose not to further complicate 
the model with these hypothetical corrections. We shall 
discuss these issues further below. 



2. The ultra-relativistic regime 



The first point to notice about eqns. H15I16|I is that 
when the contraction phase starts and the Hubble param- 
eter becomes negative the velocity will tend to increase — 
the Hubble term becomes an acceleration term, rather 
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than a damping one. This will be compounded be the 
fact that once one is beyond — 1/2 the momentmii 
term in the velocity equation also changes sign (c/. Fig. 
121), but it is clear that given the constant increase in the 
magnitude of H this term will gradually prevail (more on 
this in the following sub-section). 

Thus one can expect that, just as in the simple case of 
the circular loop, the string velocity will gradually tend 
towards unity. In this approximation, and neglecting for 
the moment the loop production and gravitational back- 
reaction terms, the evolution equation for the correlation 
length easily yields 



L (X a 



(18) 



Note that this is the same overall scaling law for the 
string network as that in a radiation-dominated expand- 
ing universe; the string network effectively behaves like 
a radiation component. In terms of the physical cor- 
relation length (|13|) we have iphys oc a, which is as if 
the strings were being conformally contracted (except for 
their rapidly growing velocities). 

However, there are several factors that must be consid- 
ered which complicate this simple state of affairs. First, 
there is the issue of loop production. Under the above 
assumptions on velocity, but putting the loop produc- 
tion term back in the correlation length equation, one 
finds the following approximate solution in the radiation- 
dominated case, 



and, in the matter era, 



(19) 



(20) 



where L„iax is the string correlation length at the time of 
maximal size of the universe, and the scale factor at that 
time was chosen to be unity (so the logarithmic correction 
in the first case is positive). Hence we see that, if c 
remains constant (or is slowly varying), asymptotically 
the scale factor dependent terms will dominate, so that 
we expect L cx In a in the radiation era, and L oc a^^'^ 
in the matter era. Again we note that the latter is also 
the scaling law for the correlation length in the matter- 
dominated, expanding universe (but again the behaviour 
of the velocities is very different in the two cases). This 
highlights the different roles played by loop production 
in the scaling behaviour of a cosmic string network in 
the radiation and matter eras, a point which has been 
noticed long ago |^ 0| in the usual expanding case. 

A strong argument can be made, however, for an im- 
portant relativistic correction to the loop production 
term in the evolution equation for the correlation length. 
In the simplest form of the VOS model there is a sim- 
ple identification between the correlation length, L, and 
the physical distance Lphys which a string segment is ex- 
pected to travel before encountering another segment of 



the same size forming a loop in the process. However, 
taking into account the Lorentz factor in H13|) , this means 
that the probability dP that a string segment will en- 
counter another segment in a time interval dt should be 
given approximately by 



dP : 



dpo 



dL 



— 1/2 

One would expect c cx 700 cx a thus driving c rapidly 
towards zero and asymptotically yielding our simple so- 
lution L (X both in the radiation and matter eras. 

Of course, during re-collapse we expect that c will de- 
pend on a number of other properties of the string net- 
work such as the enhanced build-up of small scale struc- 
ture due to the contraction. Eventually, however, the 
Hubble radius will fall below even the length scale of 
wiggles on the string after which our asymptotic solution 
L (X a? should be valid. In what follows, we shall con- 
sider the two well-motivated cases, first c — const. ^ 0, 
and secondly c = 0, the probable asymptotic case. 



,dt 



~,dt 



1/2 j- 



(21) 



3. Approaching the ultra-relativistic regime 

Returning to our analytic solutions for the constant 
loop production case (|19|l and (|20|l . we can use the ve- 
locity equation to find an approximate solution for the 
evolution of the velocity as it approaches unity. Taking 
the simplifying assumption that the momentum param- 
eter is slowly varying, one can arrive at the following 
implicit solution 



f cx a* exp 



2k{v) f a^/^da 



A 



L{a) 



(22) 



where A = 1 in the radiation era and A = 2 in the matter 
era. Substituting (|19I20|) one respectively obtains 



1 



^rad 



CX (— In a) 



ik{v) /c 



in the radiation-dominated case, and 



1 - « a'-"'"'^''/^^ 



(23) 



(24) 



in the matter-dominated case. Recall that the momen- 
tum parameter k{v) is negative in the ultra-relativistic 
regime. Hence in the limit where u — > 1 and therefore 
fc — > the asymptotic solution would have the form 



7 ^ cx (1 — v^) cx 



(25) 



The presence of the momentum corrections, which phe- 
nomenologically account for the existence of small scale 
structures on the strings, imply that convergence will be 
slower than this power. 
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4- Friction from particle scattering 

So far we have been neglecting the contribution of fric- 
tion due to particle scattering for the evolution of the 
cosmic strin g ne twork. The reason is that previous work 
[TlL Is^ im has shown that for heavy strings (say, 
those formed around the GUT scale) this will only be 
significant very close to the Hagedorn temperature, that 
is, for temperatures close to those of the string-forming 
phase transition. For the rest of the cosmic history of 
the network, this friction is always subdominant relative 
to the damping generated by the expansion itself. How- 
ever, for light strings (say, those formed around the elec- 
troweak scale) friction is significant for a longer period, 
and hence it is important to discuss how our previous 
results might be changed in this case. 

One can show 0, I3l that the effect of friction due 
to particle scattering may be characterized by a friction 
length scale 



/3T3 (3T§ ' 



(26) 



where T is the temperature (Tq ~ ImeV today), rj is the 
symmetry-breaking scale of the string and /3 is a phe- 
nomenological parameter counting the number of particle 
species interacting with the strings {(3 ^ 0{1)). Includ- 
ing this term in the string equations of motion H15I16|I is 
then fairly straightforward — see (l2Ll3i^ for a derivation. 
One finds 



2^ =2{1 + vl)HL- 



dvo 



dt 



(1 



L 



2H 



ii 



(28) 



The epoch at which the frictional force becomes subdom- 
inant is given by 



thag 



^hag 



1/2 



1 



(29) 



where M counts the total number of effectively massless 
degrees of freedom in the model (hence M = 106.75 for a 
minimal GUT model, but it can be more than an order 
of magnitude higher for particular extensions of it). 

One can then proceed as before and look for scaling 
solutions for the particular case when this extra frictional 
term dominates over the damping due to the Hubble flow. 
Again we start by considering the simplest case where 
both loop production and gravitational back-reaction are 
neglected. Assuming a relativistic network w w 1, the 
criterion for friction dominating the velocity evolution is 
> 2H which in the collapsing radiation era occurs at 
when the scale factor again reaches 1221 or alternatively 
in terms of the time before the big crunch 



(Gm) 



'"pi 



(30) 



or, still equivalently, the temperature « /m^i with 
a* « 77^/(r7ipiTo). This expectation can be confirm by 
solving the above evolution equations numerically, as 
shown in Fig. |31 For t > and while the string remains 
relativistic {v « 1), it is easy to see that exponential 
brakes are applied to the momentum with 



1 r 



jv = 7*1;* exp ; 



1 



a, 
a 



(31) 



Even if 7,?;, was extremely large at i = i*, the strings 
would soon be driven non-relativistic when the scale fac- 
tor had shrunk by only the logarithmic factor am/a* « 
l/ln(7*w*). 

From this point onwards, the velocity evolution be- 
comes analogous to friction domination in flat space, with 
Hubble anti-damping irrelevant and the small velocities 
simply reflecting a balance between the curvature scale 
and the friction length: 



L 



(32) 



When the velocity is driven away from unity w << 1 at 
tar, the previous scaling law L oa breaks down. Sub- 
stituting (|32|l into the evolution equation for the corre- 
lation length, then yields the simple asymptotic solution 
for a << ttni-, 



L cx a , 



V oc a 



(33) 



The string network ends in friction domination by com- 
ing to a standstill in comoving coordinates and then being 
genuinely conformally contracted. The energy density in 
the string network now behaves as poo oc a~^, so if it 
has not dominated the energy density before 1^1-, then 
it never will before the next bounce or final big crunch. 
This final demise for the network is expected. As the 
background temperature and density approaches those 
of the original string-forming phase transition the strings 
should effectively dissolve back into the high density ra- 
diation background from whence they had come. 



D. Further improvements 

1. Dynamical friction 

We now discuss some further contributions to the dy- 
namics of the cosmic string network in some detail, as 
well as a few caveats to this approach. An effect related 
to friction from particle scattering is dynamical friction, 
first discussed in . Cosmic string dynamics will 

be damped by dynamical friction, the magnitude of the 
effect depending both on the nature of the background 
fluid, the string mass per unit length and the amount of 
small scale structure of the cosmic string network. 

For the simple case of a wiggly string oriented along the 
z axis moving in a matter background this effect can be 
easily computed from the fact that in the rest frame of the 



8 




FIG. 3: Illustrating the effect of friction due to particle scat- 
tering in the expanding and contracting phases. A GUT-scale 
string network is evolved forward for 10 orders of magnitude 
in time starting at the epoch of formation (solid lines). The 
two lines correspond two different initial conditions, charac- 
teristic of first and second order phase transitions. At this 
point the evolution is switched to the contracting phase, and 
the network evolved back until the Hagedorn temperature, 
(dashed lines). The dotted line marks the end/beginning of 
the friction-dominated epoch, a, (see the discussion in the 
text). 



string the velocity change in the direction perpendicular 
to the string experimented by a matter particle is given 
by 0: 



5v„ 



2iTG{il - f) 



(34) 



where Vs is the string velocity with respect to the matter 
particles and fi and T are the effective string mass and 
tension per unit length, respectively. It should be obvious 
that the total momentum transferred to the particles in 
the y direction vanishes because there is a cancellation 
due to particles passing above and below the string. This 
is not the case in the x direction where we may calculate 
5vx from the fact that the value of the velocity of the 
mass particle should be the same before and after the 
interaction has taken place. This means that to second 
order in 5vy we have Sv^ = —^{6vy)'^/vs so that 



Svx = 



8Tr^G^{il-T)il 



(35) 

The total momentum transferred to a particle in the x 
direction in the rest frame of the string is then m'jSvx. 
This means that the total momentum transferred to the 
string in the interval of time dt in this frame is given 



approximately by 



SvxpVsdt 



dz 



dy ■ 



ZvsSvxdt 
327rG 



(36) 



where ^ ~ H^^ is the string length and we have as- 
sumed the universe to be flat with background density 
p — 3H'^ / {8ttG). Making a Lorentz transformation to 
the rest frame of the fluid we finally obtain that the frac- 
tion, f, of the string momentum, p = pjS.Vs, lost in one 
Hubble time is given by 



dpH- 
dt p 



37rG {jl- Tf 



4/itif 7^ 



pvs 



ip-T) 



fs7 



(37) 

We note that this simple analysis is linear on Svy which 
effectively means that we are assuming Gp <C 1 and Vs S> 
Gi~p-f). 

In the context of a collapsing universe the ultra- 
relativistic (7^1) regime is the most relevant. In this 
regime we have / ~ 0{Gp) during the matter dominated 
era. The effect of dynamical friction in a radiation back- 
ground is less dependent on the small scale structure on 
the string and is simply given bv|4(l : 



/ 



dpH- 
~di 



P 



3GA7(1 + 
327r 



(38) 



In the ultra-relativistic regime we have / ~ 0{Gpj) dur- 
ing the radiation dominated era. We note that these 
results correct equation typos (by a factor of 7) in refs. 

Hence, we conclude that in a collapsing universe and 
during the matter era dynamical friction is never able to 
slow the strings down significatively so long as Gp <C 1. 
On the other hand, during the radiation era dynamical 
friction will have a bigger impact on string dynamics and 
could in principle halt the evolution of 7 at 7 ~ {Gp)~^. 
Having said this, it is important to keep in mind that 
this analysis assumes a homogeneous and isotropic back- 
ground and so may not strictly apply in our case. More- 
over, the fact that a significant amount of the momentum 
will be transferred from the strings to the background 
(specially in the radiation era regime with 7 ~ (Gp)~^) 
will in itself add to the anisotropics which naturally occur 
in our model. 



2. A string-dominated collapsing universe 

If the energy density in the cosmic string network be- 
comes a non-negligible contribution to the overall energy 
density of the universe before tnr, then the Friedmann 
eqn. ^ must be modified to include the string den- 
sity contribution on the right-hand side. Furthermore, it 
may even happen that the string network becomes the 
dominant contribution to the dynamics of the universe, 
and in that case, the scaling laws can be significantly 



9 



modified — this type of scenario was first studied, for the 
case of expanding universes, in 

In our case, one finds that to zeroth order the 
radiation-hke behaviour of the network, given by Eqns. 
(|18I25|I . is maintained. However, if one calculates (using 
the same method as sketched above) the first order cor- 
rections to this behaviour, one finds that they are slightly 
different, namely 

L = {L„iax - clna) (39) 

and 

1 - oc a4+2M-)/A ^ (40) 

where we have defined — STrG/i/S, that is a measure 
of the string energy scale. 

Notice that even though the loop chopping efficiency c 
still affects the behaviour of the correlation length, it no 
longer affects the behaviour of the velocities (at least in 
an explicit way). In this case the corrections to the 7 cx 
scaling depend only on the momentum parameter 
k{v) and on A. 

Also notice that A is a very small number, e.g. 
/S.GUT ~ 3 X 10~^ for a GUT-scale network and /!^ew 
3 X lO"^'^ for an electroweak-scale network. The correc- 
tion term to the Lorentz factor scaling law is still becom- 
ing less and less important as the velocity increases, since 
the momentum parameter is approaching zero, but the 
fact that A is so small implies that the convergence to- 
wards the asymptotic solution is slower in this case — the 
more so the lighter the strings are. 

3. Other effects and caveats 

In the analysis in section Fill CI we have not explicitly 
included the effect of gravitational back-reaction. How- 
ever, these solutions will still hold when this is incorpo- 
rated. Indeed, in the context of the VOS model one can 
rigorously show that although this term will clearly 
affect the quantitative values of the parameters in a given 
scaling solution, as well as the timescale needed for such 
solution to be reached within a given accuracy, it can not 
affect the existence of such solutions. 

One relevant issue is that of loop reclassification. In the 
expanding phase a small loop chopped off by the string 
network will slowly decay into gravitational radiation. 
However, in the collapsing phase it is possible that an 
initially small loop becomes a large loop (that is with 
a size in the center of mass (CM) frame comparable to 
H^^) before it can lose all its energy. The size of the 
loop in the CM frame is given approximately by: 

L{t) = Lo - rG/i(i - io) 

where to is the physical time when the loop was produced, 
Lq is the initial loop size, and the numerical coefficient T 
is independently of the loop size , but does depend on its 



shape. Numerical simulations have shown that (F) ^ 65. 
The lifetime of the loop is then: 

L 

Hence, we see that for r > H^^ the loop will not have 
enough time to lose all its energy before the big crunch 
which means that loops produced with sizes greater than 
TG^H^^ will eventually re-join the long string network. 
This effectively contributes to a smaller value of c. 

A final important element not included in the above 
solution is the presence of small-scale 'wiggles' on the 
string network and the loops 0, 0| . The momentum 
parameter k{v) can to some extent account for this, but 
only in a very simple and phenomenological way. As dis- 
cussed elsewhere in this article, in the case of a contract- 
ing universe, we expect that the effect of the small-scale 
structures on the string dynamics should be proportion- 
ally larger than that in the case of an expanding universe. 
Although one could attempt to use more elaborate ana- 
lytic models to model this [l3,II2,liSlil, we shall leave 
this for further work, and for the moment restrict our- 
selves to the simple and more intuitive VOS model. 

Having said that, and even considering that the model 
has been extensively tested and shown to be accurate, 
one should keep in mind that in this case it is not ex- 
pected to do as well as in the previously studied cases. 
Apart from the issue of small-scale structures, it is also 
worth emphasizing that the VOS model assumes that the 
long string network has a Brownian distribution on large 
enough scales, which may not be a realistic approxima- 
tion in a closed, collapsing universe. This point clearly 
deserves further investigation. Indeed, given the non- 
trivial (fractal) properties of a defect network even in 
the simplest linear scaling regime, it is quite interesting 
to ascertain what are the statistical properties of the net- 
work in a closed universe around the epoch of maximal 
volume. We leave this topic for future analysis. 

E. The equation of state 

An approximate equation of state for a cosmic string 
network in the relativistic limit is easy to obtain (see for 
example Kolb and Turner [50|): 

p, = (2«,2-l)^. (41) 

This result can be obtained by taking an average over 
all possible directions of the energy-momentum tensor of 
a straight string. The generalization of this result for a 
domain wall network is also straightforward [sof: 

2 

Pw = {vi - -)p^ . (42) 

The same asymptotic limit is obtained when v —^ 1 for 
both cosmic strings and domain walls: they behave es- 
sentially as radiation. The caveat here is that this deriva- 
tion explicitly assumes a 'perfect gas' of strings or walls: 
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in other words, it assumes that there are no dissipative 
effects, either due to the defect motion or to defect inter- 
actions. Another minor caveat is that one should have 
the gas in a box much larger than the network correla- 
tion length, which will not be the case in a closed universe 
around maximum expansion, just like the Brownian as- 
sumption. 

This result could also be obtained from energy- 
momentum considerations taking into account that the 
comoving momentum should be proportional to a~^. For 
a point mass this means that jpVp ex. a~^, for a straight 
string ajsVs oc and for a planar wall a^jwVw oc a~^. 
This means that when w — > 1 we should have respectively 

7p (X (43) 
7s oc a"^ (44) 

7w oc . (45) 

Taking into account the variation of the comoving vol- 
ume in obtaining the density of each one of these objects 
during the collapsing phase it is straightforward to show 
that Pp^s.w oc which is just radiation-like behaviour. 

IV. STRING NETWORK EVOLUTION: 
NUMERICAL SIMULATIONS 




FIG. 4: The same string network box is run without (solid) 
and with (dashed) physical resolution point joining. This cru- 
cially affects the loop production. The two panels show the 
fraction of the total string energy that is in the form of long 
strings (left), and the total number of loops in the simulation 
box (right), as a function of the conformal time (defined such 
that maximal expansion occurs at 77 = 1 and the big crunch 
at 77 = 2). This clearly shows that high resolution is crucial 
in this case. 



With the caveats discussed above in mind, we are now 
ready to study the numerical evolution of a cosmic string 
network. 



A. The numerical code and basic checks 

We have performed a number of very high resolu- 
tion Goto-Nambu simulations on the COSMOS super- 
computer, using a modified version of the AUen-Shellard 
string code . 

We have simulated a simplified scenario, where a radi- 
ation or matter dominated universe is evolved in the ex- 
pansion phase for a while, to allow the initial conditions 
of the string box to be erased away (for example, the ini- 
tial string network is set up in a cubic lattice), and then 
reversed the sign of (the square root of) the Friedmann 
equation, thus forcing the universe to collapse. This set- 
ting is numerically simpler to simulate than a closed uni- 
verse, or one whose collapse is induced by a negative 
cosmological constant, but is still good enough to allow 
us to test the validity of the above solutions. 

In addition to the main set of simulations, we have 
performed a range of other control simulations to test 
various numerical issues and satisfy ourselves that our 
results had good enough numerical accuracy. Part of 
these control simulations will be briefly described below. 

We emphasize that given the expected enhancement 
of small-scale structure on the strings, having very high 



resolution is a crucial factor, especially in what concerns 
the loop population. In particular, time-saving schemes 
like point joining are not accurate enough, even if one 
keeps a constant physical resolution. This is illustrated 
in Fig. 01 which shows the outcome of evolving the same 
simulation box with and without a point joining approx- 
imation. One sees that the number of loops produced is 
very different, and an even more dramatic discrepancy 
can be noticed by looking at images of the long strings 
in the box as the simulation evolves. We believe that the 
main simulations shown below have a sufficient resolution 
(in terms of points per correlation length) and dynamic 
range to provide statistically significant results. 

The first point to be established is that there is a solu- 
tion with the rough properties of the one described above, 
namely with the velocity evolving from approximately 
V ^ l/v^ (corresponding to the usual scaling solution in 
the expanding case) to u ^ 1 as the universe collapses, 
and that such solution is stable. With this aim we have, 
among other tests, carried out some simulations which 
start collapsing after an expansion of a single time step, 
but which start out with v ^ 1. The outcome of these is 
illustrated in Fig. 

There is a very clear decrease of the velocity at early 
times, and an apparent tendency for an increase at late 
times, though it must be said that the numerical accuracy 
slightly deteriorates at the very end of the simulation. 
Notice that the decrease is much stronger for the radia- 
tion case. At the same time, the power law dependence 
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FIG. 5: Evolution of the long-string velocity (left panel) and 
scaling law of L oc a'' (right panel) for simulations starting 
with abnormally high velocity at the epoch of maximal size of 
the universe. Both are shown as a function of conformal time 
(defined such that maximal expansion occurs at = 1 and 
the big crunch at 77 = 2). Solid and dashed lines correspond 
respectively to matter and radiation universes. 

of the correlation length on the scale factor, L oc , is 
approximately constant and in the range 3/2 < /3 < 2, 
except for very early in the simulation where the net- 
work hasn't yet erased its initial conditions. Hence we 
can conclude that these results are consistent with the 
existence of an attractor solution of the type described 
above. Further evidence for this convergence has been 
observed by other means, such as starting simulations at 
maximal expansion with much lower velocity (say zero) . 

B. Detailed scaling solutions 

With some confidence in the robustness of our meth- 
ods, we can now proceed to study very high resolution 
simulations of the behaviour of a string network in a uni- 
verse with a period of expansion followed by contrac- 
tion. The result of two such simulations, for universes 
filled with radiation and matter, is shown in Fig. El and 
two snapshots (typical of the expansion and contraction 
phases) of a fraction of the simulation box in one of the 
runs is shown in Fig. \7\ 

Notice that for times up to 77 ~ 0.5 the network hasn't 
yet erased its initial (lattice) conditions, so this period 
should be disregarded when considering the scaling anal- 
ysis. We should also mention that for simplicity we are 
only plotting every twentieth time step in the simulation, 
and also we are not plotting the error bars associated with 
each data point. Depending on the quantity, these are of 
the order of 10%, and are typically larger at later than 
at earlier times. 



FIG. 6: Basic properties of cosmic string networks in the ex- 
pansion and contraction phases, as a function of conformal 
time (defined such that maximal expansion occurs at = 1 
and the big crunch at = 2). The solid represents a simula- 
tion in the matter era, while the dashed line is for a radiation 
era run. Plots successively show the behaviour of the re-scaled 
correlation length L/t (top left), the scaling law of the corre- 
lation length relative to the scale factor L on (top right), 
the velocity of the long string network (middle left), the ratio 
of the loop and long string velocities (middle right), the frac- 
tion of the string energy in the box in the form of long strings 
(bottom left) and the total number of loops (bottom right). 

During the expanding phase we confirm the usual lin- 
ear scaling regime, namely 

(X a^, f 00 = const. (46) 

in the radiation-dominated case, and 

Lf-xp (X t (X a"^/^, Woo — const. (47) 

in the matter era. 

As soon as the contraction phase starts, these laws 
are modified. As expected, the velocity starts increasing, 
and the scaling of the correlation length with the scale 
factor also drops, being approximately constant to begin 
with, and then rising slowly. One can roughly identify a 
transient scaling phase, valid in the period rj ~ 1.0 — 1.4, 
where one approximately has 

Ltrans OC a (48) 

in the radiation-dominated case, and 

Ltran. CC tt^/^ (49) 

in the matter era. (These can not be easily recovered 
by analytic methods using the evolution equations for 
the VOS model discussed above.) Unfortunately, the ex- 
tremely demanding requirements in terms of resolution of 
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FIG. 7: Two snapshots of the evolution of one of our sim- 
ulation boxes — they are typical representatives of string net- 
works in the expansion (top box) and collapsing (bottom box) 
phases, respectively. Note that these snapshots are only a 
small fraction of the total simulation box. In the first snap- 
shot the horizon is about twice as large as the box shown: in 
the second they are about the same size. 



the simulation do not currently allow us to run simula- 
tions with longer dynamic range to establish beyond rea- 
sonable doubt whether this scaling law approaches (3 = 2, 
as predicted above. However, there are strong indications 
that the networks are evolving towards this asymptotic 
regime, as shown by the relatively rapid climb of the ex- 
ponent in Fig. El 

It is clearly noticeable that the velocity rises much 
faster in the matter era than in the radiation era. It 
is also interesting to point out that during the collapse 
phase the loop and long string velocities are noticeably 
different, and this difference (which is more significant 
in the radiation than in the matter case) increases with 
time. The plot also shows an apparent difference in this 
velocity ratio in the expanding phase, but this is not 
significant. The initial lattice conditions of our simu- 
lations do tend to give equal velocities to 'long' strings 
and small loops, but as they start evolving, the fact that 
the strings are on a lattice means that there can be no 
inter-commutings for the first few time steps, and this ar- 
tificially makes the long string velocities fall behind those 
of the loops. Eventually, once the inter-commutings start 
and the network erases the 'memory' of its initial condi- 
tions, this velocity difference is also gradually erased. 

Finally, we also notice that the network keeps chop- 
ping off loops throughout the simulation, and that there 
is a dramatic increase in the small scale structure of the 
network, particularly at later times. Visually, the string 
network develops large numbers of 'knots', highly convo- 
luted strings regions where the wiggly long strings have 
collapsed inhomogeneously. These small scale features 
have proved to be difficult to evolve numerically, and this 
in fact turns out to be the main limiting factor which at 
present preventing us from running the simulations closer 
to the big crunch. 



C. Contrasting the analytical and numerical 
approaches 

We end this section by testing our solutions for the 
evolution of the string network in the contracting phase, 
obtained with the VOS model, against our numerical sim- 
ulations. The outcome of these tests is summarized in 
Fig. E 

We compare our solutions H19I2U|I for the long string 
correlation length and (|23I24|) for the long string RMS 
velocity with the numerical simulations described above, 
by plotting respectively L/a^ and (1 — v'^)/a*. We have 
assumed the usual ansatz for the momentum parameter 
k{v) — see Eqn. ((T7|l and ^3 — while for the loop chop- 
ping efficiency we have assumed a constant value c — 0.23 
which was obtained in previous studies of high resolution 
simulations in expanding universes [l3 | . 

We can see that, given the approximate nature of our 
solutions and the numerical errors in the simulations, the 
matching between the two is fairly remarkable. Naturally 
the quality of the fit will be crucially determined by our 
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FIG. 8: Comparing our analytic solutions for the behaviour 
of cosmic string networks as they approach the asymptotic 
regime with numerical simulations in the matter (solid) and 
radiation (dashed) epochs. The top panel refers to our so- 
lutions for the correlation length (I19I20I I. while the bottom 
panel refers to the solutions for the long string RMS velocity 
1231241 1. In all cases the dotted curves are obtained using the 
analytic solutions and a loop chopping efficiency c = 0.23. 



are likely to be applicable only in asymptotic regimes. 
The loop chopping efficiency is clearly important initially 
but it cannot continue to grow arbitrarily large, since the 
amount of energy transferred to loops (and small scale 
wiggles) at any time is limited by causality. On the other 
hand, a scenario in which loop production switches off 
completely may be attained ultimately during the final 
collapse. However, this is clearly delayed by the obvious 
build-up of small-scale structure on the strings. So at the 
moment, given the finite resolution of the simulations, the 
quantitative behaviour of the network is somewhat open 
to debate. Any more accurate modelling can only be 
meaningfully done in the context of a proper multi-scale 
string evolution model which, although possible 0, 0, 
li^ l is beyond the scope of the present work. 

Finally it is also worth keeping in mind that any dis- 
cussion of the evolution of a cosmic string network with 
the present formalism is only applicable while one is well 
below the Hagedorn temperature, at which the strings 
would 'dissolve' in a reverse phase transition. Discus- 
sions of asymptotic regimes should be taken with some 
caution, since a cosmic string network will only survive 
the bounce intact if this happens before the Hagedorn 
temperature is reached. 



assumptions about c and k{v), and as such one could con- 
sider this an independent test on the behaviour of these 
parameters. However, given the uncertainties discussed 
above, one can not really meaningfully use these simula- 
tions to 'measure' c directly from the simulations (with 
robust error bars) and test the k(v) ansatz. 

Among other reasons this is pertinent because, as men- 
tioned above, there could be further corrections to the 
equations in the analytic model. 

We note that tentative measurements of the loop chop- 
ping efficiency c directly from the simulations (that is, 
without using any VOS model dependent assumptions) 
are consistent with no variation throughout the range 
probed by the simulations, though again error bars in 
this particular measurement are significant. (This is be- 
cause this measurement requires the calculation of sec- 
ond derivatives of quantities in the simulation, whereas 
all quantities plotted in this paper require, at most, the 
calculation of first derivatives, and are therefore much 
more stable numerically). However, the limited dynam- 
ical range of the simulation does not allow us to distin- 
guish between the two discussed ansatze for the evolution 
of c (namely c = constant and c oc 700^^^). 

If anything, one could use the observation that in Fig. 
IHlthe numerical curves are steeper than the analytic ones 
to infer that our approximate solution is underestimat- 
ing the loop production, though on the other hand this 
could be due to further corrections on k{v) rather than 
to corrections on c. In any case, it is clear from direct 
inspection of snapshots of the simulation boxes that the 
amount of small-scale wiggles is gradually building up. 

Be that as it may, the two extreme cases outlined above 



V. DISCUSSION: COSMOLOGICAL 
CONSEQUENCES 

The evolution of the string energy density is depen- 
dent on the dynamics of the universe. In an expanding 
universe the long string energy density will evolve as 

Poo oc a" (50) 

where n « —4 during the radiation-dominated era, n « 
—3 during the matter-dominated era and n « — 2 during 
a curvature-dominated or accelerated expansion era. The 
overall density of strings remains constant relative to the 
background density p in both radiation and matter eras 

^^aGfi, (51) 
P 

with CTr w 400 and a^n ~ 60 respectively 0, 
During curvature domination or accelerated expansion, 
the string density grows relative to the other matter as 
Pool Pra OC a. For GUT-scale strings with Gp, ~ 10^^ this 
gives the interesting conclusion that today strings have 
a comparable energy density to the cosmic microwave 
background radiation Pcmb/Pm ~ 10~^. However, a real- 
istic cyclic model will continue to expand well beyond 
to, so the string density at maximum expansion will 
end up being much greater than the radiation density 
(even for strings considerably lighter than GUT scale). 
In addition, the gravitational (or other) radiation pro- 
duced through the continuous decay of the string network 
evolves as 

Pgr oc . (52) 
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It might appear that this contribution would become neg- 
ligible during the matter era but in each Hubble time the 
strings lose about half their energy into gravitational ra- 
diation, so this background always remains comparable 
to the string density pgr ~ Poo- 



A. A string-dominated universe 

Now consider the collapsing phase in which the string 
network, like the gravitational waves they have produced, 
begins to behave like radiation. Globally, the density of 
both the strings and the gravitational waves will grow 
as and, together with any other radiation compo- 
nents, they will eventually dominate over any nonrela- 
tivistic matter. In a realistic cyclic model reaching maxi- 
mum expansion in the far distant future, sufficiently mas- 
sive strings and their decay products will have a greater 
density than the cosmic microwave and neutrino back- 
grounds. As the universe contracts, then, it will eventu- 
ally reach a state in which the relativistic string network 
and/or their gravitational waves dominate the global dy- 
namics of the universe! This seems likely to lead to a 
dramatically different energy content for the universe af- 
ter it emerges from the next bounce. 

At this stage we should note, however, that if the gravi- 
tational radiation background becomes dominated by the 
longer wavelength modes then the radiation fluid approx- 
imation will eventually break down, and one then expects 
that this background will behave as 'curvature' rather 
than radiation — this is shown, in an the context of in- 
flation, in [sj- Obviously perturbations in such modes 
can't be directly detected or have a direct impact on cos- 
mological observables, but they can have an impact on 
the background in which they propagate. On the other 
hand, it is not yet clear if these results are directly ap- 
plicable to the contracting case, and if and when that 
regime is reached, since in order for the long wavelength 
modes to become dominant one requires that loop pro- 
duction (and hence gravitational radiation) switches off 
fast enough on small (sub- horizon) scales. 

Even for lighter strings which do not dominate the uni- 
verse, they would end up with a much greater density 
in the collapsing phase than they had previously during 
expansion. If the universe went through a bounce, the 
energy density in the cosmic strings and gravitational ra- 
diation produced by the network would be much greater 
after the bounce than before it (though the exact amount 
is dependent on the model details, in particular on the 
duration of the matter and curvature and/or accelerated 
expansion era). For example, bounds on the string mass 
per unit length obtained in order not to overproduce a 
gravitational radiation background may be severely mod- 
ified, in addition to the more general constraints on ad- 
ditional relativistic fluids [s^ 0| . 



B. An inhomogeneous universe 

Furthermore, unlike the uniform CMB background, 
the energy density in both cosmic strings and gravita- 
tional radiation will be very inhomogeneous. In the col- 
lapsing regime, we expect that an increasingly small frac- 
tion of Hubble regions will have a string passing through 
them. Those that do will become string dominated since 
the string energy density in those regions will approxi- 
mately evolve as 

Roc -2 fco\ 

^ cx 7oo oc a , (53) 
P 

up to the corrections described in the previous section. 
This means that for these regions the assumption of a 
FRW background will cease to be valid at sufficiently 
late times, and the defects can make the universe very 
inhomogeneous or even anisotropic Even Hub- 
ble regions without strings can be expected to have large 
fluctuations in their gravitational radiation content. For 
sufficiently massive strings, both of these effects can be 
expected to survive the bounce to create large inhomo- 
geneities in the next cycle. 

We conclude that a cosmic string network will be a 
significant source of radiation, entropy and inhomogene- 
ity which may be problematic in the cyclic context. An 
attempt to solve the problem of the overproduction of 
entropy and unwanted relics has been proposed in the 
ekpyrotic context 0- They suggest a cyclic model for 
the universe where an extended period of cosmic acceler- 
ation at low energies is used to remove the entropy, black 
holes, and other debris produced in the previous cycle. 
It is clear that when quantitative calculations are carried 
out to establish the amount of acceleration required to 
dilute unwanted debris, then the answer will depend on 
whether or not cosmic strings or other topological de- 
fects are present; a much longer period could be required 
if they are. 

We note that something analogous happens for the case 
of black holes — see for a discussion of this case. From 
cycle to cycle one expects that they will accumulate, since 
they will tend to have Hawking lifetimes longer than the 
duration of the cycle. The only ways to get rid of them 
are having a bouncing universe that is very close to flat- 
ness (so as to increase the duration of the cycle) or having 
then annihilated or torn apart at the bounce singularity. 

Finally, some of the results described in this paper are 
also expected to be valid for other topological defects, in 
particular domain walls [tI . We shall reserve a more de- 
tailed analysis of the evolution and implications of strings 
and other defects in the context of bouncing universes to 
a forthcoming publication. 

VI. CONCLUSIONS 

In this paper we have presented a first study of the ba- 
sic evolutionary properties of cosmic string networks in 
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contracting universes, using both analytic methods and 
higli- resolution numerical simulations. Wo have shown 
that the string network becomes ultra-relativistic, and 
at late times will approximately behave like a radia- 
tion fluid. We have derived new analytic scaling solu- 
tions describing this behaviour and shown, through high- 
resolution numerical simulations, that these analytic so- 
lutions are a good approximation to the actual string 
dynamics. 

The main cosmological consequence of this asymmet- 
ric behaviour in the evolution of cosmic string networks 
in the collapsing and expanding phases is that it makes 
them a significant source of entropy and inhomogeneity, 
and therefore establishes the need for a suitable entropy 
dilution mechanism if they are present in a bouncing cos- 
mological scenario. This mechanism will also operate, 
mutatis mutandis, for other stable topological defects. 
Conversely, if direct evidence is found for the presence of 
topological defects (with a given energy scale) in the early 
universe, their existence alone will impose constraints on 
the existence and characteristics of any previous phases 



of cosmological collapse. 
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